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Abstract — The Geosynchronous Imaging Fourier Transform 
Spectrometer (GIFTS) Sensor Module (SM) Engineering Demon- 
stration Unit (EDU) is a high resolution spectral imager designed 
to measure infrared (IR) radiances using a Fourier transform 
spectrometer (FTS). The GIFTS instrument employs three Focal 
Plane Arrays (FPAs), which gather measurements across the 
long-wave IR (LWIR), short/mid-wave IR (SMWIR), and visible 
spectral bands. The raw interferogram measurements are radio- 
metrically and spectrally calibrated to produce radiance spectra, 
which are further processed to obtain atmospheric profiles via re- 
trieval algorithms. This paper describes the processing algorithms 
involved in the calibration stage. The calibration procedures can 
be subdivided into three stages. In the pre-calibration stage, 
a phase correction algorithm is applied to the decimated and 
filtered complex interferogram. The resulting imaginary part 
of the spectrum contains only the noise component of the 
uncorrected spectrum. Additional random noise reduction can 
be accomplished by applying a spectral smoothing routine to the 
phase-corrected blackbody reference spectra. In the radiometric 
calibration stage, we first compute the spectral responsivity based 
on the previous results, from which, the calibrated ambient 
blackbody (ABB), hot blackbody (HBB), and scene spectra can 
be obtained. During the post-processing stage, we estimate the 
noise equivalent spectral radiance (NESR) from the calibrated 
ABB and HBB spectra. We then implement a correction scheme 
that compensates for the effect of fore-optics offsets. Finally, for 
off-axis pixels, the FPA off-axis effects correction is performed. 
To estimate the performance of the entire FPA, we developed an 
efficient method of generating pixel performance assessments. In 
addition, a random pixel selection scheme is designed based on 
the pixel performance evaluation. 

I. Introduction 

The Geosynchronous Imaging Fourier Transform Spectrom- 
eter (GIFTS) Sensor Module (SM) Engineering Demonstration 
Unit (EDU) is a high resolution spectral imager designed to 
measure infrared (IR) radiances using a Fourier transform 
spectrometer (FTS). The GIFTS instrument employs three 
Focal Plane Arrays (FPAs), which gather measurements across 
the long-wave (685-1130 cm -1 ) IR (LWIR), short/mid-wave 
(1650-2250 cm” 1 ) IR (SMWIR), and visible spectral bands. 
The raw interferogram measurements are radiometrically and 
spectrally calibrated to produce radiance spectra, which are 
further processed to obtain atmospheric profiles via retrieval 
algorithms. In this paper, we emphasize on the algorithmic 
components of the calibration process, which are divided 
in multiple stages. Prior to the radiometric calibration, the 
decimated and filtered complex interferograms must be phase- 
corrected so that their imaginary spectra contain only noise. 
The phase correction is performed on all three types of 


measurements for a series of scans, i.e., the ambient blackbody 
(ABB) reference, hot blackbody (HBB) reference, and scene 
scans. To reduce the random noise embedded in the reference 
measurements, a spectral smoothing method can be applied to 
the phase-corrected blackbody reference spectra. In the next 
stage, the calibrated ABB, HBB, and scene radiance spectra 
are computed based on the previous results. An assessment 
of the noise equivalent spectral radiance (NESR) is derived 
from the calibrated ABB and HBB references. During the 
post-processing stage, we implement two additional correction 
algorithms to compensate for the fore-optics offsets and FPA 
off-axis effects. In the last section, we introduce an efficient 
method of estimating the noise performance of the entire 
FPA, as well as a random pixel selection strategy that can 
be employed for generating a representative pixel sample set. 

II. Single-Pixel Calibration Algorithms 

For a single pixel element on the FPA, the calibration pro- 
cedures involve the following steps: phase correction, spectral 
smoothing, responsivity calculation, radiometric calibration, 
NESR assessments, fore-optics offsets adjustment, and off-axis 
effects correction. 

A. Pre-calibration Algorithms 

1) Phase Correction: The phase correction algorithm is 
applied to the decimated and filtered complex interferogram 
7 m (n). Its complex spectrum 7V m (cr) has the form of 

N m (a) = N(a)e^ +n(a), (1) 

where n(a) represents the noise vector, and N (a) is the real 
and noise-free spectrum. To estimate the phase function 0(cr), 
we adopt the Forman- Vanesse-Steel (FVS) method [1]. In this 
technique, a Hamming window w(n), of length W, is applied 
to the complex interferogram with respect to the zero path 
difference (ZPD), i.e., n = n zp d- The windowed interferogram 
is written as I w (n) = w(n)I m (n ), where 

0.54 - 0.46 cos (^zy), \n - n zpd \ < , 

0, otherwise. 

( 2 ) 

The resulting interferogram I w (n) is circular shifted to obtain 
I' w (n) such that I' w ( 0) = I w (n Z pd)- The spectral response of 
I' w (n ) is computed via the Fourier transform, i.e., A0(a) = 
Hi' w (n)}, from which, the phase function can be estimated 



LW Real (Smoothed) and Imag. Spectra Pixel (64, 65), 25 scans, ABB00494_HBB00483 .res^O.STcm -1 
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Fig. 1. Longwave phase-corrected and smoothed real and imaginary 
blackbody spectra at ambient (260K) and hot (286K) temperatures, Pixel (64, 
65), 25 scans. 


the offset term can be solved from 


N A (a)B H (a)-N H (a)B A (a) 

°( a > = =H~, - =A~ • 

N (a) — N (a) 

Next, the calibrated HBB, ABB, and scene radiance spectra 
can be determined, based on the responsivity and offset results, 
from 


- ^wr~ 0 ^ 


(9) 

( 10 ) 


respectively. The variable i denotes the scan number, i.e., i = 


1 , 2 ,. 


,S. 


as (j>(a ) = tan 1 By applying the phase estimate, 

the phase-corrected spectrum becomes 

N(a) = + n(a)e ~ j ^ a) , (3) 


and its resulting real and imaginary components are given as 
Re{7V(cr)} ~ 7 V(<t) + Re{n(cr)} cos(0(<j)), (4) 

Im{7V(cr )} « Im{n(cr)} sin(0(cr)). (5) 

Note that the imaginary part of the spectrum contains only 
noise, which can be discarded. 

2) Spectral Smoothing: The phase correction is performed 
on three types of spectra, i.e., ABB, HBB, and scene. An 
optional smoothing algorithm can be applied to the blackbody 
references to improve the noise performance of the calibrated 
radiances. Two smoothing techniques are being considered 
for this purpose: the least-squares smoothing filter and the 
moving average (MA) filter. The MA smoothing method 
approximates the function within the predefined window as a 
constant whereas the least- squares algorithm fits a polynomial 
to all points within the window. The MA method does not 
preserve higher order moments, however, it is less demanding 
computationally. See [2], [3] for detailed reviews of these 
two techniques. When the final noise performance of these 
two methods are compared, no significant differences were 
observed; therefore, the MA smoothing technique is chosen 
due to its simplicity. Fig. 1 illustrates the phase corrected 
and smoothed longwave spectra for the real and imaginary 
blackbody references at ambient and hot temperatures. 


B. Radiometric Calibration 

To obtain the calibrated radiances, the first step involves the 
calculation of the spectral responsivity, which can be computed 
from 


R(a) = 


—H . 


N H (a) - N A (a) 
B H (a) - B A (o) ' 

-=A, 


( 6 ) 


(<t) = |Ef=i *?(*) and N (a) = | Ef=l 

represent the mean HBB and ABB spectra averaged over S 
scans, respectively. Using the relation of 


tH 


N (a) = R(a)(B H (a) + 0(a)), 


C. Post-processing 

1) Noise Equivalent Spectral Radiance (NESR): This sec- 
tion describes the estimation of the noise equivalent spectral 
radiance (NESR) from calibrated blackbody references. The 
NESR is generally considered as a measure of the instrument 
noise performance, and can be estimated as the standard devi- 
ation of calibrated blackbody radiance spectra from multiple 
scans. The NESR definition is given by 


NESR H ’ A (a) 


\ 


IZWm 




B H ’\a)) 2 . 


( 11 ) 


Fig. 2 depicts the smoothed HBB and ABB NESR estimates 
over their corresponding non-smoothed NESRs computed 
from 25 interferogram scans. 


LW ABB and HBB NESR Estimates Pixel (64, 65), 25 scans, ABB00494_HBB00483 



Fig. 2. Smoothed and non-smoothed longwave HBB and ABB NESR 
estimates computed from 25 interferogram scans. 


2) Fore-optics Offsets: The previous calibration results are 
computed using internal blackbody references. To compensate 
for the offsets generated from fore-optics, a correction scheme 
is incorporated into the previous model. This can be achieved 
by using additional data measurements collected from an 

H a 

extended source. Given C e (a) and C e (a) as calibrated 
extended blackbody sources averaged over all scans, we have 

Cf’ A (cr) = R e (a)B^’ A (a) + O e (a). 


(7) 


( 12 ) 


By solving for R e (a) and O e (cr), the actual scene radiances 
can be approximated as 




Nfja) - O e (<7) 

-Re (O') 


(13) 


Fig. 3 shows the calibrated longwave HBB, ABB, and scene 
radiance spectra after performing the fore-optics offsets cor- 
rection. 
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spectra for two pixels, one pixel is located near the center 
of the FPA whereas the other is situated near the comer of 
the FPA. A better wavenumber alignment was achieved after 
performing the off-axis effects correction. 


LW Calibrated Spectra (Off-Axis Corrected) for Pixels (64,65) and (1,1), Scan 1 , CCE00006 .res^O^cm" 1 



Fig. 4. Calibrated spectra for Pixel (1,1) and Pixel (64, 65) after performing 
the off-axis effects correction. 


Fig. 3. Calibrated longwave HBB, ABB, and scene radiance spectra after 
the fore-optics offsets correction (25 scans). 


III. Focal Plane Array (FPA) Performance 
Assessment 


3) Off-axis Effects: In this section, the correction of the 
FPA off-axis effects is presented. Since the GIFTS FPA 
contains 128 x 128 detector elements, the distortion caused by 
the off-axis effects is too great to be neglected. The off-axis 
pixels are sampled at slightly shorter optical path differences 
(OPDs), which cause the spectra of these pixels to expand 
to slightly higher wavenumbers [4]. The correction can be 
formulated as a fractional sampling rate conversion problem, 
which can be solved via sine interpolations, i.e., 

M—l 

N i (<r) = ^2 ( ct/ ) sinc ( a - a ')> a ' = 7 ’ ( 14 ) 

where a is the on-axis wavenumber scale, a' is the off-axis 
scale, and / is the off-axis factor that is computed from the 
FPA geometry. However, the sinc interpolation calculation is 
computational intensive, furthermore, the truncation window 
associated with the sinc kernels may cause overshoots in 
the resulting data. It has been shown that the zero-padding 
interpolation in the interferogram domain is equivalent to 
the sinc interpolation in the spectral domain except that the 
zero-padding method cannot evaluate the interpolation output 
at an arbitrary point [5]. To overcome this limitation, we 
designed the method of “over-padding”, in which, we assign 
an over-padding factor g that is closely correlated to the final 
resolution of the corrected spectra. For instance, if / = 0.9977, 
and the length of the spectrum is N, then the zero-padded 
interferogram length is the roundoff value of g * N/f. If 
g = 100, the actual off-axis factor of 0.997702 is obtained; 
if g = 1, then the actual value of / becomes 0.99758 due 
to roundoff errors. The corrected spectrum can be obtained 
by downsampling the Fourier transform of the over-padded 
interferogram by the factor g. Fig. 4 depicts the corrected 


We have discussed the algorithmic procedures required for a 
single pixel measurement. To evaluate the entire FPA without 
completing all of the calculation steps, we designed a pixel 
inventory method, which estimates the FPA responsivity and 
noise distributions. 


A. Pixel Inventory 

The responsivity for any detector element can be approxi- 
mated from the normalized intensity value at the ZPD, which 
is given by 


R x ,y = 


I (n Z pd) 


j^EE^'W 

x y 


V x,y, 


(15) 


where x , y are pixel numbers. The FPA responsivity estimate 
is shown in Fig. 5(a). The FPA responsivity distribution within 
a predefined range of 0.8 < R x,y < 1.2 is plotted in Fig. 5(b). 
To assess the FPA noise, each interferogram is first normalized 
by its value at n = n zp d , i.e., 

- T x,y (n) 

I x J{n)= ™ 1 L Vs,y. (16) 

J-m \n zl pd) 

Next, the estimated noise is computed from the root-mean- 
square value of the last L samples in each interferogram that 
is represented by 


jyx,y = 


\ 


1 M—l 

l \lm y (n)\ 2 V x,y. 


(17) 


Similarly, the FPA noise and its distribution are shown in 
Figs. 6(a) and 6(b), respectively. 





Fig. 5. Longwave pixel inventory: (a) FPA responsivity estimates, and (b) Fig. 6. Longwave pixel inventory: (a) FPA noise estimates, and (b) FPA 
FPA responsivity distribution plot. noise distribution plot. 


B. Random Pixel Set Generation 

Once we have the knowledge of the FPA responsivity and 
noise estimates, a set of selection criteria is established by 
defining an acceptable range for the responsivity and noise. 
Thus, a random sample set of pixels can be chosen such that 
it reflects the statistical properties of the entire FPA. We can 
then perform the actual calibration procedures on this sample 
set. Fig. 7 shows a sample longwave pixel selection scheme, 
in which 64 random pixels are selected with 4 pixels per tap. 
The selection is uniformly distributed and unique within each 
tap. 


IV. Conclusion 

This paper describes the required stages for calibrating 
the GIFTS SM EDU data spectrally and radiometrically. The 
following procedures were introduced: phase correction, spec- 
tral smoothing, radiometric calibration, NESR computation, 
correction algorithms for the fore-optics offsets and the off- 
axis effects. In the last section, methods for estimating the 
FPA’s performance and generating random pixel sample sets 
were discussed. 
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